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Abstract The current structural models of the cellulose 
microfibril as well as its mechanical and thermal properties 
are reviewed. The cellulose microfibril can be considered as 
a single thin and long crystalline entity with highly aniso¬ 
tropic physical properties. The contribution and limit of 
different methods employed such as electron microscopy, 
infrared spectroscopy, X-ray scattering and diffraction, 
solid state nuclear magnetic resonance spectroscopy, and 
molecular modeling are also discussed. 
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Introduction 

The fine structure of native cellulose fibers has important 
implications in the biological processes of plant growth, 
chemical or enzymatic transformation of cellulosic biomass, 
and for the conception of new material and understanding 
of physical properties of a wide range of cellulose-based 
materials. A considerable part of cellulose research has 
been devoted to its elucidation. Diverse models and struc¬ 
tural concepts have been developed based on different 
observations and scientific environments, some of them 
facing counter evidence but still persisting in the literature. 
Even the term “microfibril” can sometimes generate confu¬ 
sion, because different authors use the term in different 
ways with closely related terms such as “elementary fibril,” 
“protofibril,” or “nanofibers.” Although the object under 
scrutiny has nanometric lateral dimensions, the term “cel¬ 
lulose microfibril” contains notions that are not limited to 
its size, as will be seen in the following. For this reason this 
historical term is maintained throughout this review. The 
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common features and variability of the basic component in 
native cellulose, the microfibril, are outlined with basic 
observations in support of the model, as well as the limit of 
the present knowledge. 


Morphology 

The cross section 

Oriented native cellulose fiber samples such as ramie give 
strong X-ray diffraction spots in the direction perpendicular 
to the fiber axis. The three strongest diffraction spots cor- 

o 

respond to spacings of about 3.9, 6.0, and 5.4 A (Fig. 1), 

o 

where the 3.9-A spacing gives the most intense spot. This 
indicates that the cellulose chains are straight ribbons 

o 

aligned in the planes, spaced at 3.9 A. Because all the 
hydroxyl groups are equatorial to the pyranose ring plane, 
it is surmised that hydrogen bonding occurs in this plane; 
this is referred to as the “hydrogen-bonded plane” in the 
following discussion. Although the approximate molecular 
conformation and the approximate lateral size of the crys¬ 
talline region and their orientation have been resolved 
rather directly by using X-ray diffraction, the organization 
of the crystallites inside a fibrillar morphology observed 
with the electron microscope has remained hypothetical for 
a long time. 

In the unicellular giant green alga, Valonia , sharp X-ray 
diffraction patterns can be obtained. The orientation of the 

o 

crystallographic plane with the 6.0-A spacing lying parallel 
to the cell wall was reported by Sponsler. 1 Although this is 
not the plane where hydrogen bonding is expected, the 
tendency to orient this plane parallel to the surface has been 
observed in hydrolyzed wood cellulose, 23 and thus a near- 
rectangular section with the crystallographic surface paral- 

o 

lei to the plane with 6.0-A interplane spacing was conceived 
as the building block. 4 The interchain distance inside this 

o 

plane is about 5.4 A, and is the shortest among different 
crystallographic planes; it is referred to here as the “dense 
plane.” Based on molecular mechanics calculations, the 
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Fig. 1 . Bright-field diffraction contrast image of a section of Valonia 
ventricosa (left, courtesy of Dr. Chanzy) and a molecular model of 
cellulose I (right) with the monoclinic (black line) and triclinic (thick 
gray line) unit cell 


formation of this plane was considered to be the most 
energetically favorable especially in a high-dielectric con¬ 
stant medium 5 such as water, and thus to precede the crys¬ 
tallization event. However, microfibrils having better 

o 

perfection in the 5.4-A-spaced “second dense plane,” can 
also be found in different genera in the order of Zygnema- 

A r-j Q 

tales, such as Spyrogira , Micrasterias , or Closterium. In the 
cell walls of higher plants such as birch 9 or flax, 10 the orien¬ 
tation of the cellulose crystallites is random in a plane per¬ 
pendicular to the fiber axis. 

The imaging of a single crystalline domain by using a dif¬ 
fracted electron beam gave the first direct proof that the 
crystalline domain extended over the whole width of the 
microfibril. Furthermore, crystalline lattice has been 

1 r\ IQ 

directly observed for Valonia , Cladophora microfibrils, 
and hydrolyzed fragments of ramie fiber, 14 extending over 
the whole width of the observed microfibril. Lattice images 
in the cross section of Valonia cell wall, oriented 
microfibrillar bundle in glomerulocytes of tunicate, ' or acid- 
hydrolyzed tunicin 17 give directly the cross-sectional shape 
of the crystals, because the two-dimensional lattice can be 
obtained only when the chain axis is almost parallel to the 
electron beam. The Valonia cellulose has a square-shaped 
cross section with the “dense plane” and the “second dense 
plane” exposed (Fig. 1), whereas tunicate cellulose has a 
parallelepipedic shape with “dense plane” exposed as the 
largest facet. The shape also coincided with the earlier dif¬ 
fraction contrast observation and the envelope of microfibril 
observed by negative staining or heavy metal shadowing. 
Thus, the microfibrils in these highly crystalline samples 
essentially correspond to the single crystalline entity. 

The concept of microfibrils with a single crystalline core 
has been further strengthened by considerations of the bio¬ 
synthesis mechanism. Transmembrane globular structures 
arranged in two dimensions, called terminal complexes 
(TCs) have been observed at the fracture surface of the 
external plasma membrane associated with the impression 
of microfibrils in different organisms. The arrangement of 


TCs had close correlation with the morphology of microfibrils 
produced by the organism, and thus the TCs would simul¬ 
taneously synthesize many cellulose chains that crystallize 
into a microfibril. In the case of vascular plants, TCs are 
considered to have a six-membered rosette morphology on 
the internal side of the membrane. These rosettes have 
been observed in a number of organisms associated with 
cellulose production. In the case of mung bean root, immu- 
nogold labeling proved that the rosettes had a high similar- 
ity with the cellulose synthase in cotton, and thus the 
rosette TC is widely accepted as the cellulose synthase 
complex in higher plants. 

Due to the small size and the tendency to tightly associ¬ 
ate with each other in the cell wall, the morphology of the 
microfibrils in higher plants is less clear. The cross-sectional 
morphology of the microfibril is considered to be close to 
square or parallelepipedic based on electron micrographs 
of negatively stained sections of populus tension wood 20 
and electron diffraction contrast images of a cross section 
of flax fiber 10 where square dots can be observed. The 
lateral dimension of the dots corresponded roughly to the 
crystallite size estimated from X-ray diffraction line 
broadening. 

The estimation of the cross-sectional dimensions with 
higher precision is problematic. Atomic force microscopy 
(AFM) is considered to have high vertical resolution, and 
thus is potentially capable of precise measurement provided 
that microfibrils can be sufficiently dispersed on a flat 
surface. This aspect of AFM seems to be underexploited 
despite the widespread availability of the technique. Many 
other techniques have been employed to evaluate microfibril 
width but most are based on nontrivial assumptions. 

Line broadening of X-ray diffraction (XD) 

The relation between crystal size and the line broadening 
of the X-ray diffraction, provided that the broadening due 
to the instrument is negligible, is straightforward. The peak 
profile should just follow the Fourier transform of periodic 
structure with limited size. However, in practice, the mea¬ 
surement of peak width is not obvious, except in a special 
case where the crystallite size is large enough compared 
with the unit cell dimensions so that peaks do not overlap. 

o 

Such cases occur for Valonia or tunicin, and the 3.9-A plane 
of ramie cellulose, which were the first to be measured by 
this method. 21 

For other cases, profile analysis based on more or less 
arbitrary intervention of the experimentalist is required, as 
was illustrated by Hindeleh and Johnson. 22 They used a 
third-order polynomial background function and a mixture 
of Lorentzian and Gaussian functions as the peak function 
to fit the equatorial intensity profiles. Depending on the 
mixing ratio of Lorentzian and Gaussian functions, the 
lateral width of the crystallites varied from 3.35 nm, when 
integral breadth was taken for 100% Lorentzian peak, to 
6.7 nm when peak width at half maximum of a Gaussian 
peak was taken. The integral breadth should give a weight 
average of crystallite size provided that the background is 
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properly subtracted. It should be noted that there is nothing 
to justify the third-order polynomial or to prefer any other 
background function. The peak function also depends on 
the crystal morphology, their size distribution, as well as 
internal strain and imperfection. Thus, the measure is highly 
model dependent, and one has to be cautious when compar¬ 
ing the values among different authors. 

An important discrepancy exists, even for well-defined 
thick microfibrils such as those from Valonia , between 
microscopic observation and lateral crystallite size esti¬ 
mated from line broadening. The standard value from 
microscopic observation is twice as large as the one deduced 
from X-ray analysis. 2 ' 

Small-angle scattering and negative-staining 
electron microscopy 

Small-angle X-ray scattering (SAXS) was first applied by 
Heyn 24 to estimate the lateral size of the crystallites in highly 
oriented fibers. Among different swelling conditions, 
straight Guinier plots, expected from rod-like morphologies 
free of interparticle correlations, were obtained for jute 
swollen in 5% sodium hydroxide or ramie and flax swollen 
in 1% alkali in the q range of about 0.015-0.1 A 1 . The 
diameter of crystallites, assimilated to cylindrical objects, 

o 

was estimated from the slope of the plot: for jute, 28 A; flax, 

o o o 

28 A; ramie, 43 A; cotton, 55 A. The range of diameter cor¬ 
responded to the results obtained from electron micro¬ 
graphs of sections of the same samples negatively stained 
using uranyl acetate. 25 However, interpretation of electron 
micrographs must proceed with caution due to the overlap¬ 
ping of layers, 26 or due to phase-contrast artifacts. 22 

A similar study was performed on Picea abies applying 
SAXS 27 combined with negative-staining transmission elec¬ 
tron microscopy (TEM). 28 In this case, the scattering profile 
was fitted in the range 0.15-0.45 A~ . The minimum of inten¬ 
sity observed at 0.3 A 1 was interpreted as a sign of rather 
narrow distribution in diameter, in contrast to the wider 
distribution obtained by TEM ascribed to the grain size of 
the contrasting agent. However, it has to be noted that the 
intensity rise after the minimum coincides with the tail of 
the first diffraction spot at 1 A 1 , and thus to interpret the 
apparent minima as a consequence of homogeneous cross 
section can be misleading. 

Microbeam SAXS of a flax fiber could be fitted by a 

o o 

mean diameter of 19 A with a standard deviation of 9 A, or 
a volume average of 30 A, in the 0.04-0.4 A 1 range. 29 These 

o 

values are similar to the 28 A measured by Heyn, but are 
significantly smaller than the crystallite width estimated 
from the line broadening of XD, as mentioned by the 
authors, even when taking the lowest estimation. Although 
the homogeneous cylinder model can explain the SAXS 
feature, a unique solution cannot be derived if the system 
is polydisperse and departs from the cylindrical morphol¬ 
ogy. The fact that the estimated diameter was smaller than 
the smaller limit of crystallite size estimated from XD raises 
serious doubts on the validity of the method in estimating 
the microfibril width. 


Quantification of surface chain 

Other methods are based on the hypothesis that all acces¬ 
sible or amorphous features correspond to the surface 
chains. The downfield shift of C4 in solid state nuclear 
magnetic resonance (NMR) spectra at 90 ppm is character¬ 
istic of crystalline cellulose, compared with 80 ppm for cel¬ 
lulose in solution. When the crystallinity of cellulose is 
destroyed by ball milling, the peak occurs at around 85 ppm. 
A similar situation can be seen for C6 where native cellu¬ 
lose gives a chemical shift of 66 ppm, while ball-milled cel¬ 
lulose has a peak at 63 ppm.' 0 

By investigating cellulose from different sources, Earl 
and VanderHart 31 noticed a correlation between the pro¬ 
portion of surface chain to chain inside the crystal, as esti¬ 
mated from the microscopically determined microfibril 
size. Wetting the sample resulted in resolution of the 
higher magnetic field components, which exhibited at least 
two distinct peaks, as can be seen in the spectra recorded 
by VanderHart and Atalla. 32 Newman 33 focused on this 
NMR component, having well-resolved peaks at 84.5 and 
83.5 ppm in never-dried kraft pulp, and ascribed it to the 
surface of the crystallites. 

Wickholm et al. 34 obtained the “surface” signal for highly 
crystalline Cladophora cellulose at 84.5 and 83.5 ppm (Fig. 
2). Because the latter resonance was more intense, they 
tentatively assigned the 84.5 and 83.5 ppm signals to the two 
exposed crystallographic surfaces corresponding to the 
“second dense plane” and the “dense plane,” respectively. 34 
In support of this assignment can be found the following 
observations: 

1. A similar difference in intensity is observed on bacterial 

cellulose (Fig. 2b). 32 



Fig. 2. Solid state 13 C nuclear magnetic resonance spectra of C4 region 
of wet native cellulose, a , Cladophora ; b , bacterial; c, cotton. 
Spectra b and c were shifted upheld by 1 ppm to match with more 
recent spectra 
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2. The 83.5 ppm component is better defined than the 84.5 
ppm component in the tunicate, 34 which does not expose 
the “second dense plane.” 16,17 

3. The peak at 84.5 ppm is more pronounced with samples 
obtained by glycerol annealing of highly crystalline cel¬ 
lulose III where the “second dense plane” is visibly 
better developed in the diffraction profile. 

4. The longer relaxation time of the 84.5 ppm component 
(18 s) compared with the 83.5 ppm component (11 s) is 
also in agreement with molecular dynamics calculation 
that predicts longer relaxation time for the “dense plane” 
surface.' 6 

The densely packed sheets to behave in a more rigid 
fashion is also informative. However, further interpreta¬ 
tions of the data require various assumptions on the peak 
profile and would encounter a similar problem to that dis¬ 
cussed for XD analysis. Even in “highly crystalline” cellu¬ 
lose, the microfibril cross section displays rounded corners. 
Thus, a surface constituted purely of two crystallographic 
planes seems to be an oversimplification. 

In the same manner, the state of surface chains can be 
estimated by the hydrogen-deuterium exchange rate of 
hydroxyl groups using either infrared spectroscopy 37,38 or 
neutron inelastic scattering measurements. 39 Both methods 
give rough agreement with the generally accepted microfibril 
size, but further ambiguity remains on the evaluation of the 
depth of the accessible hydroxyl group at a given condition. 
For infrared investigation, the difference in molar absorp¬ 
tion of the OD/OH stretching bands for the crystalline core 
and surface is another nontrivial parameter. 

Longitudinal direction 

In highly crystalline cellulose from algae or tunicate, and in 
bacterial cellulose, it is relatively easy to obtain microfibrils 
sparsely dispersed to follow individual microfibrils over a 
long distance. Dark-field or bright-field diffraction contrast 
images often show patches along the microfibril. It needs to 
be considered that the diffraction contrast depends on the 
Bragg condition and thus these images are sensitive to the 
angle of the crystalline plane with respect to the incident 
electron beam. These patches often move continuously 
along the microfibril during the observation probably due 
to the inclination drift of the supporting membrane. Thus, 
the microfibril is of single crystalline nature over a long 
distance along the microfibril, even if twists or strain are 
also omnipresent. 

In higher plants, one of the first microfibrils to be 
observed in the dispersed state, were those of quince slime, 
in which glucuronoxylan surrounds the surface preventing 
the microfibrils from sticking together. 40 This allowed visu¬ 
alization of an individual microfibril over a long distance, 
sometimes exceeding 6 pm. The majority of the microfibril 
had a width in the range of 3-5 nm, but some of them were 
thinner than 2 nm. 41 

Secondary-wall cellulose such as that of wood cellulose 
needs some chemical treatment to allow good dispersion. 
Lepoutre and Robertson 42 succeeded in obtaining well- 


dispersed microfibril of wood pulp by heterogeneous 
modification of never-dried wood pulp and measured the 
microfibril width to be in the range of 3.5-4.5 nm. More 
recently, selective oxidation of primary hydroxyl groups on 
the surface of microfibrils using never-dried samples, 43,44 or 
carboxymethylation of never-dried pulp solvent-exchanged 
in isopropanol 45 were used to facilitate the separation of 
microfibrils. In all these cases, thin continuous microfibrils 
with constant lateral size, at least within the resolution of 
TEM and AFM, of a few nanometers (<5 nm) over a length 
of a few microns can be seen with occasional kinks. Thus, 
as a basic rule, microfibrils can be considered as a crystalline 
entity continuous over a very long distance, having a range 
of lateral dimension depending on the biological origin and 
development stage. 46 

Periodic defects along the microfibril have been expected 
from the acid hydrolysis behavior of cotton fibers. Although 
“amorphous regions” and “surface chains” on the microfibril 
share many spectral features,' 9 the measurement of chemi¬ 
cal reactivity shows the rupture of hydrogen bonding 
between 03 and the ring oxygen only in this region of 
defect. 47 Recently, small-angle neutron scattering revealed 
a Bragg diffraction peak for partially deuterated ramie 
fiber, indicating a periodical distribution of regions suscep¬ 
tible to deuterium substitution. The weak contrast seen in 
neutron scattering experiments and the small weight loss 
during acid hydrolysis indicates that this susceptible part is 
extremely small compared with the periodicity of 150 nm 
and represents a few residues every 300 residues. 48 

The origin of the periodical defects along the microfibrils 
in higher plants is not clear. An interesting observation is 
reported on cellulose microfibrils produced in vitro from 
membrane fragments, which cannot be fragmented by acid 
hydrolysis, whereas the in vivo counterpart will be frag¬ 
mented as any other cellulose from higher plants. 49 The high 
density of microfibrils in the cell wall would result in a tight 
lateral aggregation' concentrating the internal strain in a 
limited zone distributed along the chain axis. The strain 
accumulates to such a level to induce even the rupture of 
the 03H-05 hydrogen bonding. 47 


Crystal structure of cellulose \ a , \ p , and IV, 

As mentioned in the previous section, all native celluloses 
share a common feature when looking at strong diffraction 
at low resolution. Small differences had been noticed 
between bacterial cellulose or Valonia cellulose, initially 
named type A, and higher plant cellulose, named type B, 
when diffraction spots with smaller intensities or infrared 
spectra were analyzed. 50 The diffraction from Valonia cel¬ 
lulose, 51,52 although very sharp and informative, could only 
be indexed with a large unit cell comprising eight chains. 
Solid-state L 'C NMR analysis showed that all cellulose 
spectra could be explained as a linear combination of two 
subspectra with different combination ratios. Diffraction 
analysis from a selected area of microfibril containing the 
pure crystalline form then allowed establishment of two 
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simple unit cells: a triclinic and one-chain unit cell and a 
monoclinic two-chain unit cell. 54 Each of them contains two 
crystallographically independent glucosyl residues. These 
unit cells imply that the basic difference between the two 
structures is how the consecutive hydrogen-bonded planes 
lie on top of each other. If the crystallization is preceded by 
a lateral association in the dense or second dense plane, the 
l a and Ip is determined already at this first event. 

Some native celluloses, often giving only blurred diffrac¬ 
tion due to the small crystallite size, are assigned to another 
unit cell called cellulose IV! similar to 1^ but with the mono¬ 
clinic angle close to a right angle. Due to the peak overlap 
and broadness, the determination of peak position is quite 
delicate. Furthermore, the diffraction pattern is the convo¬ 
lution of a Laue function that comes from the lattice repeat, 
and the molecular form factor inside the lattice. As a con¬ 
sequence, the experimental peak positions are influenced 
by the slope of the form factor. 55 Cellulose in red meranti 
has a rectangular unit cell, taking the crystallite size effect 
into consideration. 56 The cellulose IV! is very close to cel¬ 
lulose Ip in many features, and the experimental protocol to 
convert cellulose III to cellulose IV : generates cellulose 1^ 
of small crystallite size when a highly crystalline sample was 
employed. 35 However, whether the “real” unit cell in 
samples classified as cellulose IV : systematically deviates 
from the cellulose 1^ needs more careful examination. 

Parallel-up structure 

Although the approximate orientation of the pyranose ring 
inside the unit cell is rather trivial from the intensity ratios 
of the three most intense reflections, whether the chains are 
packed in a parallel or antiparallel manner is less clear. 
X-Ray fiber crystallography assisted primitive molecular 
modeling approaches in the 1970s and favored a parallel 
arrangement. 52,57 59 However, the structures proposed inde¬ 
pendently by two groups at the time had opposite chain 
direction inside the unit cell when a common convention 
was used. French et al. 60 compared different refinement 
methods and datasets employed by different authors, and 
concluded that the main limit in discriminating “up” and 
“down” structure lied in the accuracy of the dataset rather 
than the computational program utilized. The difference 
between the “up” and “down” structure is essentially the 
different relative displacement of the hydrogen-bonded 
planes. A direct indication for a parallel structure was then 
obtained by selective staining of the reducing end and 
directional enzymatic hydrolysis. 62 These methods, however, 
did not indicate whether the structure was parallel-up or 
parallel-down. This ambiguity was lifted by combining the 
reducing end labeling and electron crystallography on a 
single fragment: the chain polarity could be deduced to be 

S'! 

parallel-up according to the standard convention. 

The diffraction data could be further improved by using 
uniaxially oriented specimens composed of relatively large 
crystallites found in nature. The diffraction data extended 

o 

to “atomic resolution” of 1 A with about 300 measurable 
intensities. 64,65 The direct localization of electron density in 


reasonable position using a partial model for phasing indi¬ 
cates a high reliability of the method. The parallel-up struc¬ 
ture was confirmed in this study. 

Hydrogen bonding 

Direct evidence of the existence of hydrogen bonding can 
be readily seen in infrared spectra, which have been pub¬ 
lished since the 1950s with special emphasis on hydrogen 
bonding. 50,66 An empirical relation exists between the hydro¬ 
gen bond distance and the position, width, and intensity of 
OH stretching bands that shift to lower wavenumber when 
an acceptor is close. 67 This red shift of the OH stretching 
band is qualitatively explained as elongation of the OH 
bond. However, because the cellulose crystal has six inde¬ 
pendent OH groups in a unit cell, which can couple with 
each other to give in-phase and out-of-phase bands, the 
assignment of the spectra can be complex. A general remark 
can be made that the cellulose I gives stronger parallel 
bands. This indicates that native cellulose has relatively 
strong hydrogen bonding along the chain axis and relatively 
low intermolecular polar interaction. This aspect is quite 
different from cellulose II and cellulose III, which exhibit 
OH bands with important red shift perpendicular to the 
chain direction and small red shift parallel to the chain 
direction, indicating a more extensive intermolecular hydro¬ 
gen bonding. 

Further insight or resolution can be obtained by external 
stimulation such as stress or temperature. Hinterstoisser et 
al. 68 70 used a two-dimensional correlation representation 
using stress variation on cellulose sheets and discussed on 
molecular mechanism of load bearing with respect to the 
hydrogen bonding. 71 However, the measurement was done 
on only slightly oriented sheets formed by a dynamic former, 
and thus the stress conditions on each microfibril are far 
from clear as will be seen in the following section. Highly 
resolved OH stretching can also be obtained by highly ori¬ 
ented thin films constituted of relatively thick cellulose 
microcrystals. 72 Marechal and Chanzy 7 ' studied highly crys¬ 
talline cellulose I p, presenting difference spectra with tem¬ 
perature modulation and the effect of heavy water wetting 
and drying. It is noted that absorbance builds up at around 
3445 cnT at the expense of the absorption at 3240 cm - , 
which disappeared on heating or decreased on drying. This 
led to the assignment of the most intense and parallel band 
at 3340 cm 1 to the 03H secondary alcohol hydrogen- 
bonded to the ring oxygen, because it should be quite stable. 
The band at 3240 cnT 1 was assigned to a minor population 
of 02H that strongly hydrogen-bonded to 06 of the subse¬ 
quent residue. This hydrogen bonding was considered to be 
easily ruptured by heating or drying. 

The location of hydrogen atoms related to hydrogen 
bonding is difficult to determine by X-ray diffraction tech¬ 
niques: the hydrogen atoms contribute very little to the X- 
ray scattering, and the hydrogen of the hydroxyl groups are 
labile. Neutron diffraction has a great advantage in that 
neutrons interact with hydrogen or deuterium at similar 
amplitude as carbon and oxygen, and hydrogen and deute- 
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rium contribute with inverse phase. Thus, deuterium sub¬ 
stitution of hydroxyl groups gives a substantial modification 
in neutron diffraction intensity. The hydrogen position 
clearly showed up for 03H in the Fourier difference maps 
using the neutron intensities from deuterated and hydroge¬ 
nated samples. However, other hydroxyl hydrogens were 
less clearly localized on the Fourier map, suggesting a 
hydrogen-bonding disorder. 64,65 The existence of two coex¬ 
isting hydrogen bonding systems, one major and one minor, 
was proposed based on structure refinement against neutron 
data. Both forms of hydrogen bonding lie on the “hydro¬ 
gen-bonded plane” where the major one corresponds to 
what was present in most structural models since the scale 
model. 74 The minor component was refined to a hydrogen 
bonding having similar feature for I a and 1/3, but it cannot 
be excluded that many different hydrogen-bonding schemes 
exist in reality, as was proposed based on molecular model¬ 
ing. 75 The structure did not significantly change at low tern- 
perature, ' because the hydrogen bonding did not converge 

r l r l 

to a single hydrogen-bonding system even at 15 K. 


Elastic properties 

The linear elastic behavior of a material can be approxi¬ 
mated by a second-order tensor, relating the stress tensor 
and strain tensor. In the system with lowest symmetry, the 
triclinic system, the number of independent coefficients is 
21, whereas in the monoclinic system with P2 : symmetry, 
the number is 13. However, as will be described in the fol¬ 
lowing, experimental values are sparse, partly due to practi¬ 
cal difficulties related to the small crystal size. 

The Young’s modulus of the crystalline domain has been 
measured by applying uniaxial tensile load along the fiber 
axis and monitoring the lattice deformation along the chain 
axis using X-ray diffraction. This is supposed to correspond 
to the inverse of a diagonal element of the compliance 
tensor. Most of the studies have been done on ramie fibers 
due to its low microfibrillar angle and single fiber length 
that allows easy manipulation. The stress on each crystallite 
is in general considered to be identical to the stress applied 
to the macroscopic sample. The values based on this assump¬ 
tion fall in the range of 115-140 GPa. 78 82 The argument in 
support of this “series model” (uniform stress) is that the 
crystalline modulus measured at 9% and 32% moisture 
content was 130 GPa whereas the macroscopic modulus was 
24 and 12 GPa, respectively. 79 However, the uniform stress 
condition is somewhat contradictory to the current model 
of the microfibril where long crystalline domains exist along 
the chain direction with very little, if any, amorphous 
domain (Fig. 3). If a parallel model is taken, the estimated 
crystalline modulus would be higher 82 because the cross 
section of crystal that bears the load will be smaller. 

In a series of Chamaecypans obtusa (Hinoki), the ratio 
of [crystal strain]/[surface strain] was reported to be highly 
correlated with the microfibrillar angle, and consequently 

QQ 

with the tensile modulus along the fiber direction. For 
wood specimens with small microfibrillar angle, the surface 



Fig. 3. Series model (left) and microfibril model (right) of a fiber under 
tensile stress. The amorphous region or surface chains (gray lines) is 
expected to have lower elastic modulus than the crystalline core (black 
lines). In the series model, all stress is transmitted through crystalline 
regions, whereas in microfibrils the crystalline core and surface chains 
share the stress 


strain was very close to the macroscopic strain. The stress 
condition would be easier to understand with smaller 
microfibrillar angle, in which case the deformation clearly 
follows a parallel situation (uniform deformation). The 
origin of high macroscopic strain of ramie cellulose com¬ 
pared with crystalline strain remains unclear, although 
analysis based on combination of parallel and series models 
can be made. 80,82 

Another experimental value for tunicate cellulose, 143 
GPa, 84 is based on Raman band shift. Strain was applied on 
two-dimensionally random oriented sample embedded in 
epoxy resin, where the stress was estimated from the shift 
of the 1095 cm" band. The empirical relationship between 
the band shift with applied stress, common to native and 
regenerated cellulose 85 was used in this case. The estimation 
assumes crystal strain to be equal to macroscopic strain and 
that the efficiency factor that relates Young’s modulus of 
the individual fiber to the modulus of two-dimensional 
sheet applies at this scale. 

The lattice deformation can also be measured in the 
lateral direction during uniaxial tensile stretching. This is 
related to the off-diagonal element of the compliance 
tensor. The strain normal to the hydrogen-bonded plane £ a 
can be measured relatively easily because the correspond¬ 
ing diffraction peak is strong and isolated. The ratio v ac = 
eje c where e c is the strain along the chain direction, does 
not depend on the assumption of stress condition or mea- 
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surement of cross-sectional area. The reported values are, 
however, associated with relatively high estimated errors 
(standard deviation): 0.56 ± 0.1 86 or 0.38 ± 0.04 87 for ramie 
and 0.46 ± 0.1 86 for flax. 

In the case of wood pulp, a negative apparent v ac was 
measured after kraft cooking, 88 whereas intact wood 
specimens always showed positive values. 89 Due to the 
microfibrillar angle, the internal stress conditions are not 
sufficiently evident to speculate and needs further under¬ 
standing of basic properties of the microfibril. 

The theoretical elastic modulus can be calculated with 
small computation resources, provided that the force field 
and the structure are realistic. Unfortunately, up to now 
there is no force field able to reproduce the unit cell dimen¬ 
sions of cellulose allomorphs to experimental precision. It 
is claimed that the mechanical properties can be properly 
estimated even with a noncomplete force field, 90 however, 
the calculated elastic properties rely heavily on the starting 
structure model. 91 Even the contribution of hydrogen 
bonding to the stiffness varied among studies from 60% 
to less than 20%. Thus, the quantitative understanding of 
molecular behavior related to elastic deformation is still a 
long way away. 

A direct measurement of the Young’s modulus can also 
be made by using the relation between sound velocity and 
density. The X-ray or neutron scattering can probe the 
density-density correlation in time and space, from which 
sound velocity in the lattice can be extracted. The signal is 
dominated by the direction with strong structure factor, and 
thus the Young’s modulus in the direction perpendicular to 
the hydrogen-bonded plane of 14.8 ± 0.8 GPa was calcu¬ 
lated with relatively high precision. In the chain direction, 
the phonon dispersion could only be measured close to the 
diffraction position, leading to lower precision, but a value 
of 220 ± 50 GPa, significantly higher than any other experi¬ 
mental value, is reported. 9 ’ 


Thermal behavior 

The measurement of anisotropic thermal expansion is more 
straightforward because it is a first-order tensor, and the 
temperature is uniform in general. Thermal expansion coef¬ 
ficients have been measured for chain direction on hemlock- 
spruce and lateral direction on hemp from room temperature 
to 200°C, 94 lateral direction on cotton, 95 on tunicin 96 and 
Cladophora , 97 and in both directions for tension wood of 
Populus maximoviczii . 98,99 Thermal contraction was observed 
for chain direction in hemlock-spruce, but the tension wood 
showed a positive, though small, thermal expansion. The 
thermal expansion is the greatest in the direction perpen¬ 
dicular to the hydrogen-bonded sheets, of the order of 100 
ppm/K when the average is taken between room tempera¬ 
ture and 200°C. Between -150°C and room temperature, 
the value is 70-80 ppm/K for cotton 95 and 46 ppm/K for 

r-t S' 

tunicin. The tension wood showed a negative thermal 
expansion along the b -axis (in the hydrogen-bonded plane) 
whereas tunicin showed a positive value, suggesting a size 


effect on the thermal property of microfibrils. The thermal 
expansion anisotropy is well reproduced in molecular 
dynamics simulations. 100 The highly anisotropic elastic prop¬ 
erty and low thermal expansion coefficient along the chain 
direction result in a low planar thermal expansion when 
microfibrils form a two-dimensional network. 101 
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